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Abstract 



We consider two models for a pair of interacting particles in a random 
potential: (i) two particles with a Hubbard interaction in arbitrary dimen- 
sions and (ii) a strongly bound pair in one dimension. Establishing suitable 
correpondences we demonstrate that both cases can be described in terms fa- 
miliar from theories of noninteracting particles. In particular, these two cases 
are shown to be controlled by a single scaling variable, namely the pair con- 
ductance gi- For an attractive or repulsive Hubbard interaction and starting 
from a certain effective Hamiltonian we derive a supersymmetric nonlinear a 
model. Its action turns out to be closely related to the one found by Efetov 
h-^. ■ for noninteracting electrons in disordered metals. This enables us to describe 

0^ , the diffusive motion of the particle pair on scales exceeding the one-particle 

jo localization length L\ and to discuss the corresponding level statistics. For 

tightly bound pairs in one dimension, on the other hand, we follow early work 



by Dorokhov and exploit the analogy with the transfer matrix approach to 
quasi-ld conductors. Extending our study to M particles we obtain a M- 
particle localization length scaling like the Mth power of the one-particle 
localization length. 
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I. INTRODUCTION 

The concept of one parameter scaling pioneered by Thouless [TJ and generalized by 
Abrahams et al 0,[| has been very fruitful in the understanding of the localization problem 
for noninteracting electrons in a random potential. The relevant scaling parameter is the 
size-dependent conductance g(L) (in units of 2e 2 /h) determined by the so-called /3-function 
(3(g) ||. This /3-function which contains all the relevant information has in principle to be 
calculated from a microscopic theory. However, the general field theoretical formulation in 
terms of a non linear a model j||| shows the universal behavior of (3(g) and that it only 
depends on the symmetry class and the space dimension. In 2 + e dimensions the /3-function 
is accessible by diagrammatic perturbation theory f^]], and in quasi-one dimension even 
"exact" solutions are known, implying exponential localization 0,[J. Furthermore, exact 
expressions for the first two moments of the conductance have been derived |],[nj . 

In one or quasi-one dimension there is, apart from the a model, another powerful theory 



1 1] to describe the scaling behavior of the conductance. In this approach the distribution of 
transmission eigenvalues is determined by a Fokker-Planck equation, the Dorokhov-Mello- 
Pereyra-Kumar (DMPK) equation, describing how the transmission properties evolve with 
the length of the wire. In the unitary symmetry class the entire distribution of transmis- 
sion eigenvalues [|T^j and the first two conductance moments |]H| are known for arbitrary 
length scales. Moreover, it has become clear recently that the DMPK equation and the a 
model formulation of quasi-ld transport are mathematically equivalent [P~5|-|HS|]. In both 
approaches, one-parameter scaling emerges in a natural way. 

Not long ago, Shepelyansky |16| studied the localization problem of two interacting parti- 



cles (TIP). He considered a one-dimensional disordered potential and a local (Hubbard-like) 
repulsive or attractive interaction U. Assuming some generic properties of the interaction in- 
duced coupling matrix elements between localized one-electron product states, Shepelyansky 
found that certain TIP-states are extended on a scale given by the two-particle localization 
length L2 oc L\, which can be much larger than the one-particle localization length L\. Us- 
ing similar assumptions, Imry |17| was able to confirm Shepelyansky 's original result for L2 



by using the Thouless |IJ block scaling picture. He introduced the notion of a two-particle 
conductance g2{L) and assumed that it obeys the same scaling behavior as the one-particle 
conductance. In subsequent work, the enhancement of the TIP localization length was con- 



firmed by numerical studies for the transfer matrix [IS], the TIP-Hamiltonian [19] and the 



TIP-Green's function [20 



In this paper we show that both the a model and the Fokker-Planck approach can be 
succesfully applied to the TIP-localization problem, thus relating the TIP-problem to an 
equivalent problem without interactions. This allows us to identify a certain pair conduc- 
tance #2 associated with TIP-diffusion and TIP-localization and to show that gi is governed 
by the same scaling function as the one-electron conductance. Furthermore, we relate #2 to 
the statistical properties of the TIP energy levels and to the transmission eigenvalues of the 
corresponding scattering problem. 

In section [TI| we consider, generalizing Shepelyansky, the TIP-problem in d dimensions 
with a Hubbard-like repulsive or attractive interaction U. We first motivate an effective 
random matrix model that incorporates precisely the 2(i-dimensional lattice structure of 
the two-electron pair states. The only assumptions concern the statistical independence and 



the typical size of the interaction-induced coupling matrix elements between these states. 
The next step is to map the random matrix model onto a supersymmetric non linear a 
model. The corresponding action is after appropriate identifications formally identical with 
the action of Efetov's "standard" a model || for disordered conductors. This allows us to 
determine in arbitrary dimensions the Breit-Wigner form for the local density of states, the 
different frequency regimes for the diffusion of interaction-assisted pairs at scale L > L\, 
and the two point correlation for the corresponding TIP-levels. We relate g<i to a pair 
Thouless energy E^ defined as the energy range for finding Wigner-Dyson level statistics. 
However, we also obtain a subtle modification with respect to the one-particle problem: the 
diffusion is logarithmically suppressed and the size of the interaction-assisted pairs increases 
logarithmically with time. These results are in agreement with recent findings of Borgonovi 
and Shepelyansky [pH] , p^1 based on numerical and qualitative arguments. A summary of our 
main results has already been published in a separate publication [^] and we present here 
a comprehensive derivation of the different steps. 

In section |H], we consider a different TIP problem, which was first studied by Dorokhov 



24[ : a tightly bound pair of particles in a strictly Id random potential. In analogy to quasi 
one-dimensional conductors, we treat this problem in terms of the DMPK equation. The 
scattering channels for the transfer matrix formulation are now provided by the internal 
quantum modes of the bound pair. The conductance gi is then related to the transmission 
eigenvalues corresponding to these scattering channels. To identify #2 i n this approach, 
we have to determine one microscopic length scale, namely the effective mean free path or 
equivalently the localization length as a function of the microscopic parameters. Again we 
find that the "conductance" #2 for tightly bound pairs has the same scaling with system size 
as the conductance of non interacting electrons in quasi one-dimensional wires. Extending 
this approach to M bound particles in one dimension, we obtain a derealization effect which 
increases with M as L^ 4 if the particles are forced to stay within a distance of order L\ from 
each other. The derealization of certain M particle states by a local interaction will be 
discussed in a greater generality elsewhere |25| . 



II. TWO PARTICLES WITH LOCAL ATTRACTIVE OR REPULSIVE 
INTERACTIONS IN D-DIMENSIONS 

We consider two interacting electrons on a lattice in d dimensions with site diagonal 
disorder and nearest neighbor hopping elements. We assume that all one electron eigenstates 
are localized with a typical localization length L\ ^> 1 (all length scales will be measured 
in units of the lattice constant). This requires the disorder to be sufficiently weak in one or 
two dimensions, while the disorder strength should be slightly above the critical disorder in 
three dimensions. The localized one electron states <p p {r) are labeled by an index p describing 
the "center" of the state. We consider the case of symmetric two-electron wave functions 
?P( r i, r 2) = VK r 2) r i) ( e -g- both electrons have opposite spin) and express the Hamilton 
operator in the basis of symmetrized products of one-electron states \p1p2 >, given by the 
wave function [w( r iVp 2 ( r 2) H-^ Pl (r 2 )^p 2 (ri)]/\/2: 

H = 12 ( £ Pi + £ P2)\PIP2 >< P1P2I + 12 QpiP2 P3P 4 \P1P2 >< PZPi\ , (2.1) 

Pl,P2 Pi, P2, PS, Pi 



where Q PlP2 p 3P4 are the interaction matrix elements. The quantities e p are the one-electron 
energies of the localized states <p p (r). For a local Hubbard interaction with interaction 
strength U 

H int = Uj2\rr><rr\ (2.2) 

r 

we have for instance 

QpiP2P 3 Pi=< Pip2\H int \p 3 p 4 >=2Uj2 l f P i(r)v P2 (r)!f P3 {r)(f P4 {r) . (2.3) 

r 

These interaction matrix elements become exponentially small if two of the "positions" pj 
are far away on a scale much larger than the one-electron localization length. Shepelyansky 
divided [HJ (in one dimension) the problem described by the Hamiltonian ( |2.1| ) into two 
blocks: a diagonal one corresponding to matrix elements < PiP2\H\p3P4 >, where at least 
two of the positions pj are separated by more than L\, and a random band matrix with 
a superimposed strong diagonal (SBRM), where every pj is at less than L\ from the three 
others. The ensuing mathematical random matrix model has been investigated in |2T)|-p 



The two crucial simplifications for this picture are first to neglect all exponentially small 
interaction matrix elements Q plP2P3Pi ( off-diagonal terms of the first block, off-diagonal terms 
outside the band of the second block, coupling terms between the two blocks) and second to 
assume that the remaining coupling elements (|2.3| ) are independent random variables with 
typically the same size. Using the ansatz 

vp{r) ~ c^W e " |r_pl/Ll a » ' ( 2 - 4 ) 



and applying a central limit theorem argument in (|2.3| ), one may estimate [ 16 , T/T | the coupling 
strength between well coupled product states of the SBRM-block to be 

'<&»**) ~ U2L i M ■ ( 2 - 5 ) 

In (|2.4j ), we assume for simplicity that a p (r) is a random variable with (a p (r)) = and 
(a p (r)a~ p {r)) = S pp S r f. 

A. Effective Random Matrix Hamiltonian 

In this study, we will consider an effective random matrix Hamiltonian (ERMH), which 
relies only on the second of the above-mentioned simplifications. In other words, all product 
states \p1p2 > will be taken into account, and we keep the total Hamiltonian, and not only 
the SBRM sub-block where the states are strongly affected by the interaction. However, the 
coupling matrix elements (|2.3|) of the total Hamiltonian matrix will still be assumed to be 
independent (Gaussian) variables. This latter point is indeed quite a serious simplification 
since the correlations between different coupling matrix elements will be neglected. In par- 
ticular, the exact operator identity Hf nt = UH int [see fl2.2p ] is violated. However, we believe 
that our ERMH will nevertheless be quite realistic due to the strongly fluctuating diagonal 
elements, which tend to eliminate the effect of the correlations. 



We consider a 2<i-dimensional configuration space (lattice) with two coordinate vectors 
R = p\ + p 2 and j = (pi — p2)/2 corresponding to twice the center of mass and half the 
distance of the two electrons. This choice ensures that the components of these coordinates 
differ for adjacent lattice sites just by +1 or —1 if the corresponding other coordinate is kept 
constant. The ERMH consists of two parts 

H = f) + (, (2.6) 

where the 

ifi) R i = V R S RA 6 n (2.7) 



correspond to the diagonal entries e Pl + e P2 of (|2.1| ). For simplicity, we choose for the 
rf R independent random variables with the distribution function p§{r\). We have typically 
po{rj) ~ l/(2Wb) for \r]\ < W\, where Wb is twice the bandwidth of the disordered one-electron 
Hamiltonian. The operator £ corresponds to the interaction matrix elements ( p2.3[ ) . The 
matrix elements are independent Gaussian random variables with zero mean and variance 

(Kl) 2 ) = ^ (1 + ^^ )a(|jR "^ l)t;(j>(J) • (2 ' 8) 

The smooth functions a(\R — R\) and v(j) decay exponentially on the scale L\. We assume 
the typical behavior: 

a(\R-R\) ~ { : T9 Zs d -,^ur, ' 7 2 ~7 ' (2-9) 

and 

vU)-{ "_,,,.,, r.;: 7 (2.10) 



U 2 Lf d , \R- 

U2 L fd e -2\B-R\/ Ll ^ R _ 


-R\<L, 
-R\^>Lx 


f 1 , IjI^^i, 





which is justified by Eqs. ( |2.3| ) and (p.4[). The function a(\R—R\) describes how the coupling 
strength decreases exponentially as the distance of the centers of mass increases, while v(j) 
describes how the "size" of the product states reduces the coupling. 

In Eqs. ( p.9| ), (|2.10| ), we have applied the simplified estimate (|2.5| ) for the typical strength 
of the interaction coupling matrix elements. However, the presented random matrix model 
and its subsequent treatment is given in terms of the general functions a(R) and v (j). It 
is straightforward to use different, more "realistic", estimates instead of (|2.5| ), e.g. to take 
effects of ballistic or diffusive motion for length scales L < L\ into account. 

B. Nonlinear a Model 

Our goal is to investigate the spectral statistics, the transport and the localization prop- 
erties of the ERMH ( [2.6| ) applying the supersymmetry technique |5,2!J. In the following, 



we choose a particular realization of the diagonal elements rf R and perform the ensemble 



average only with respect to (. This particular average is denoted by (• • -)^. We consider 
the generating functional 



F(J) = If Dip exp 



l$(E-H + {^ + ie)A + j)il> 



= Udet' 1/2 ( y E-H+{- + ie)A + j)\ . (2.11) 

Here, if) is a supervector with components if>j (R) , which are 8-dimensional supervectors with 
entries zi, z%, Xi, Xi, z 2, %2, X2, X2, where z u (%„) are complex bosonic (fermionic) variables. 
The diagonal matrix A = diag(l4, — 14) describes the grading into advanced and retarded 
Green's functions, u is a frequency, J is a source matrix, and if) is given by if) = ip^A. For 
the graded determinant ("sdet") and for the graded trace ("str") we use the convention 
that str(A) = tr (Abb) — tr(A FF ) and sdet (A) = exp(str \n(A)), where A BB (A FF ) is the 
boson-boson (fermion-fermion) block of the supermatrix A. 

Using suitable choices |29|j30|] for J and taking derivatives of F(J) with respect to J at 
J = 0, one can obtain moments of the Green's function. These contain all the informations 
about the spectral statistics and transport properties we are interested in. Our strategy 
is the following: We derive a nonlinear a model from the functional F(J) and discuss its 
physical implications by comparing our result with the a model description of noninteracting 
electrons derived by Efetov ||. 

The ensemble average with respect to ( is readily performed and results in 

F(J) = f Dif) exp -ip(E-f) + (-+ie)A + j)ip--^a(\R-R\)stT[K(R)K(R)] 
J 2 ^ 2 ' 8 

R,R 

(2.12) 
where K(R) are 8x8 supermatrices given by 

K{R)=Y,v{j)^{R)i) 3 {R) ■ (2-13) 



We see that the statistical properties ( |2.8|) of £ allow for a convenient decoupling into the 
center of mass coordinates R and the relative coordinates j. The quartic term in fl2.12| ) 



is decoupled in a standard way [^,^9j by a Hubbard-Stratonovich transformation. This 
introduces an integral over a field of 8 x 8 supermatrices cr(R) with the same symmetries as 
K(R), and we get 



F{J) = f Daexp(-d[a}) 



^W} = ^Y,( A ~ 1 )RR st MR)^(R)} + \^UR) ln(E-f) + (%+ie)A + J+±a) . (2.14) 

-rt,ix 

Here, A -1 is the inverse of the matrix A R ^ = a(\R — R\) and a is given by 

(&) R i = v(j)a(R)S RR 8 fj . (2.15) 



The supertrace ("str^)") appearing in ([2.1 4|) extends over the full 8L 2d -dimensional super- 



vector space, where L is the system size assumed to be much larger than L\. To proceed, 
we apply as usual a saddle point approximation, which is valid in the limit L\ ^> 1. The 
corresponding saddle point equations 8C\/8cr{R) = can be written in the form (for J = 

and uj = 0) 

a(i?) = -£a(|i?-i?|)$>(j)- X . (2.16) 

r i E-r 1 3 R + ieA + \v{j)a{R) 

The sums extend over L\ d (contributing) terms and, applying a central limit theorem ar- 
gument, we may replace the particular values of rf^ by an average with respect to the 



IR 



distribution po( 7 ?) ( see below for a more precise justification). Then ( |2.16|) can be solved by 



a homogeneous (i.e. ^-independent) ansatz a(R) = T +iTiA, where T , 1^ are determined 
by the implicit equation 

r„ + iTiA = -.Bb 5>0') [ dri poirj) „ . . , r rvr , r M , (2-17) 

j J E -r] + ieA+ ±v(j)(T + iT 1 A) 

B = J2a(\R\) . (2.18) 

R 

In the limit L\ 3> 1, To and Y\ are much smaller than the typical integration range 2W& for 
r\ and we find 

To^-BoSxV fd V po(v)^— , (2.19) 

Ft-irBoStpoiE) , (2.20) 

with S u = J2j[ v {j)] u ~ Ef and V J drj (• • •) a principal value integral. From Eqs. ( |2.9|) and 
( |2.1U| ) we find the estimate 

r x ~^Li d <z:w b (2.2i) 

w b 



if U and W b are of the same order of magnitude. The imaginary part of the r.h.s. of ( 2.16|) 
is a sum of Lorentzians with a typical width Fiv(j). The most important contributions in 
this sum arise from \j\ < L\ and \R — R\ ^$ L\ such that v(j) ~ 1. The maxima rf R of the 
Lorentzians have thus a typical distance ~ 2WbLi 2 <^i T\ and we are in the regime of well 
overlapping resonances. Therefore it is indeed allowed to replace ( |2.16j ) by its average with 
respect to ?7~. As usual, the saddle point T + iT x A is not unique and we have to consider a 
complete saddle point manifold T Q + iTiQ with Q = T _1 AT, T being an element of a certain 
supermatrix group f5,29]. The supermatrix Q fulfills the nonlinear constraint Q 2 = 1 and 



is in our case a member of the orthogonal a model space. In the following, we replace <j(R) 
by To + iTiQ(R) where each Q(R) fulfills this constraint but may slowly vary in i?-space. 
The matrix A RR = a(\R — R\) has eigenvalues a(k) = J2R e * kRa (\R\) an d the eigenvalues of 
A^ 1 are just given by 



oB 2 r .. I 



W)^ B ° + k m ' ^-m?^™ • (2 ' 22) 



Here we have applied an expansion with respect to small momenta k since mainly the slow 
modes will give the relevant contributions f5,31]. In analogy to 1 31], we evaluate the first term 



of the action ( |2.14j ) by performing a Fourier transformation, applying ( |2.22| ) and afterwards 



transforming back to position space. Since Q(R) varies slowly with R, we may go over to a 
continuum limit and obtain the following expression for the functional F(J): 

F(J) = J DQ exp(-£ 2 [Q]) , 

C 2 [Q] = -j^ J dR str ([V B Q(R)] 2 ) + istr 07?) ln(E - fj + (f + ie)K + J + \a) , 

&)& = vV)F* + KiQ(R)\&RR*i3 ■ ( 2 - 23 ) 



For convenience we have kept the discrete notation for R in the second term of the action 
€2- Expression ( |2.23| ) defines the nonlinear a model suitable to describe the ERMH. 



C. Breit-Wigner Form for the Local Density of States 

As a first application, we determine the local density of states, 

pUE) = — lim Im(< R]\(E - H + ie)- x \R] >) . (2.24) 

This quantity can be calculated using small 4x4 supermatrices cr(R) in Q2.14J ) and replacing 
A by 1. Then only the saddle point r +iFi (and not a full saddle point manifold) contributes. 
Choosing for J a supermatrix that has its only nonvanishing diagonal entry xL g at site (J, R) 
(L g is a diagonal supermatrix with entries +1 and —1 in the boson-boson and fermion- 
fermion blocks, respectively), one finds 

fS R (E)= 1 - ,■ 1 ' riU0 ' ) ; ■ (2.25) 

This is a Breit-Wigner function with the width ~F\v{j). The same result can also be 
obtained by a standard diagrammatic expansion in powers of (. The complex self-energy is 
then given by by S^ = — |(r + iTi) v(J) with T and r x determined by ( |2.17| ) and fl2.18| ). 



To interpret (|2.25|) we note that the noninteracting product states with eigenenergies 
rf R acquire a finite width Tiv(j) due to the interaction. This width is of the order Ti for 
\j\ < L\ and decreases esponentially as Y\t~^' Lx for \j\ ^> L\. This means that those 
states rf R constructed from well separated one-electron states (|j| ^> L\) constitute good 
"quasi-levels" with exponentially large life times. We may consider a critical pair size L c 
such that the pair states with 2\j\ ^ L c in a system of size L SJS are still reasonably well 
coupled. For this, the effective level spacing A e g = Wb (L c L SJS )~ d has to be compared with 
the smallest possible level width Yi e~ 2Lc / Ll . This yields the rough estimate L c ~ Li\nL sys . 
The r] R with \j \ > L c contribute to a discrete point spectrum whereas the states with \j \ < L c 
are well mixed and form effective pairs of maximal size L c ~ L\ lnL sys . 



D. TIP-conductance g 2 at scale L > L\ 



We now want to describe the dynamics of these pairs in terms of the effective a model 
( p.23| ) with the action ^[Q]- In order to understand the relation with a disordered metal, 
we recall a few results for the corresponding o model description given by Efetov || and 
characterized by the action 

£ met [Q} = ~vV J dr sir ((VQ) 2 )+^-(-i^ + e) J drstr(Q(r)A) . (2.26) 

Here T> denotes the classical diffusion constant and v the local density of states. Efetov 
derived || this action for frequencies \u\ <C 1/r where r is the elastic scattering time 
due to the disorder. The action C met [Q] therefore describes the quantum diffusion for time 
scales larger than r. The nontrivial quantum properties are due to the nonlinear constraint 
Q 2 = 1 of the supermatrix Q(r). One may parameterize || Q by Q = Aj^-p where P 
is a supermatrix with only offdiagonal entries in the A-grading, APA = —P. Expanding 
the action (|2.26|) with respect to P, Efetov showed that a perturbative evaluation of the a 
model is completely equivalent to the diagrammatic theory of disordered metals in terms of 
(interacting) Diffuson and Cooperon modes. The leading order corresponds to the classical 
diffusion and the higher order terms yield the first quantum corrections. Furthermore, Efetov 
found that the action ( |2.26 ) leads to the well known level statistics [32 of random matrix 



theory if the considered range of energies is smaller than the Thouless energy E c = V/L (i.e. 
\u>\ <C E c ). For frequencies u such that E c <^ \u}\ <^1/t and in the metallic regime E c 3> A 
(A denotes the level spacing), the corrections to the universal Wigner-Dyson statistics have 
been calculated by Altshuler and Shklovskii 



|33|j . In one dimension, the a model ( |2.26| ) 
describes a (thick) wire. All the wave functions are localized [f^H with a localization length 
f = 2vrz/D (for L > £)■ 



Comparing the first terms of the actions of (|2.23|) and (|2.26|) , one gets for two interacting 
particles 



u e ffV, 



cff 



1 T\B 2 



TV 



R 2 



7rp (E) 2 SfB 2 



El 



T 2 



^2.27) 



where the second identity results from Eq. ( |2.20| ). We have introduced an effective (pair) dif- 
fusion constant V e s and an effective local density of states z/ e fj. For (quasi) one-dimensional 
systems the quantity in ( |2.27| ) is proportional to the localization length L 2 for the center of 

(U 2 /W 2 )L\. 



U 



mass coordinate for the pairs and we recover Shepelyansky's result 
We reiterate that we have taken all, also the weakly coupled, product states into account. 
It has been suggested |26| that these weakly coupled states lead to a logarithmic correction, 
L 2 ~ L\j\yyL\. From our analytical considerations we find no evidence for such a behavior. 
In order to extract the diffusion constant X> e g from ( |2.27| ), we still need a reasonable 
estimate for the effective local density of (product) states v eS that contribute to the diffusion. 
Naively, one could consider the effective size L c ~ L\ In L sys of the pair states, which we 



estimated in the above discussion subsequent to Eq. ( F2.25|) . This would yield the estimates 



iy eS ^ Po (E)L d c ^ Po (E)L d 1 {\nL 



sysj 



V, 



elf 



u 2 

W b l K 



sysj 



(2.28) 



However, this simple picture assumes that all product states up to a size of L c contribute 
equally to the diffusion. Obviously, this is not the case. Similarly to the frequency condition 
\cj\ < 1/r for a disordered metal, we have to require |u;| < l/r c to justify ( |2.28| ), where 
l/r c = Txv(L c ) = A cS is the effective level spacing of the well coupled product states, 
which is by definition equal to the inverse life time of pairs of size L c . We recall that 
l/r c = W b L~ 2d ~ W b (Li \i\L sys )~ 2d is much smaller than the energy scale I\ ~ (U 2 /W b )Li d , 
which is the inverse life time of the best coupled states (of size \j\ ^ L\ — L c /lnL sys ). For 
frequencies between l/r c and r 1; we thus expect that the number N(u) of contributing 
states depends on to via: 

N(uj) = [L cS (uj)} d , L eS (uj) = ^\n(^\ , (2.29) 

where L eS (uj) is the pair size such that the corresponding inverse life time Ti v(L eS (uj) J 
equals \u>\. Eq. ( |2.29| ) results in a frequency dependent density of states and diffusion 
constant, 

i/effM~Po(£) (^i ln ( r i/l^l)) =* ^cffM~— Lf- d (in^VJlul))' . (2.30) 

This estimate is valid for l/r c < |o>| < Ti, whereas for \u>\ < l/r c Eq. ( |2.28| ) applies. The 
estimates (|2.28| ) and ( |2.30|) , based on the simple but physical life time argument, can also be 



obtained in a more formal way directly from the a model action (|2.23|) . For this we rewrite 
C 2 [Q] (at vanishing source term J = 0) as 

C 2 [Q] = -^oiiPeff / dR str ((VQ) 2 ) + JdR f R (Q(R)) , (2.31) 

where 

f R {Q) = \ E str H E - ^ + (! + ^) A + i v U) ( r o + rtiQ)) ( 2 - 32 ) 

is some effective Q-potential for the a model. In Appendix [A|, this expression for the effective 
potential is simplified to give 

f R (Q) ~ -i^cu h (^j p Q {E) str(QA) , (2.33) 

where h{T\/u) is defined in Eq. ([A(]) of Appendix |A] and has the limiting behavior described 
in ( |A7| - |A9|) . Comparing (|2.33| ) with the second term in the action (|2.26| ), we obtain a complex 
effective density of states, 

u eS {u) = p (E)h(^j . (2.34) 

For the two cases ( |A8|) and ( |A~9|) we directly recover the estimates ( f2.30| ) and ( |2.28| ), respec- 
tively. On the other hand, the frequency range |u;| ^> Ti correponds to 

/ fl (Q)~ Jristr(QA), (2.35) 

10 



which means that the action C2[Q] no longer depends on u>. 

To summarize, we have found that for frequencies smaller than r 1; the a model action 
£2[Q] can be written in a similar form as the action £ me t[Q] for a disordered metal, 

Gift?] = -| ^(u) V cS {uo) J dr str ((VQ) 2 ) + ?^M(-i| + £ ) J dr str(Q(r)A) . (2.36) 

It was necessary to introduce both a frequency dependent diffusion constant V e s(uj) and a 
frequency dependent density of states u c s(u), 

£ b LJ- d (\nL sys y d , M<1/t c 
ZWM ~ { (2.37) 



n 6 



L\- d (Mgr))" , i/r c <|^|«r 1 



g eff r 2 £? 2 _^_ 2 



^ e ff(o;) = — , . , o- cff = — -5- ~ tj-2 -^i • (2.38) 



Here, we have defined in analogy to a disordered metal via the Einstein relation a "pair con- 
ductivity" a c g which does not depend on frequency. This determines the TIP- conductance 
gi mentioned in the introduction through the usual relation between conductance and con- 
ductivity: g% = o~ e gL d ~ 2 . The analogy with a disordered metal allows us to draw some 
very interesting and far reaching conclusions. First, we can identify the coupling constant 
|z/ eff £> eff = !<7 e fT in C2[Q] as a universal scaling parameter. The corresponding scaling func- 
tion is precisely the same as that of a disordered metal provided the latter is described by the 
"standard" a model Q2.26|) . In particular, the perturbative evaluation of the /3-function in 



2 + e dimensions 0-0] is also valid for our problem of diffusing or localized electron pairs (the 
second term in £2 ft?] containing the supertrace str(QA) is not affected under the renormal- 
ization ||). In one dimension, we are able to directly identify the pair localization length 
L>2 ~ o"cff ~ (U 2 /W 2 )L 2 . These remarks rigorously justify Imry's application [HJ of the 
Thouless scaling picture |TJ. 

E. TIP-Dynamics for L% « L < L 2 

A second conclusion concerns the pair-dynamics in the metallic (or diffusive) regime 
(L < L 2 for d — 1,2) where the a model can be treated perturbatively (cp. Ref. f|). The 
relevant diffusion propagator TZ(q,u>) = [cr eS q 2 — iuv e ff(u)] _1 = {u cfi (u)[V cS (u!) q 2 — iuj]}^ 1 
now contains the frequency dependent diffusion constant ([2.371) an d density of states Q2.38|) . 
This gives rise to some subtle modifications of standard diffusion. The precise calculation 
of the diffusion propagator lZ(R,t) = (2n)~( d+1 ^ f dcu f dq 1Z(q,uj)e i<yqR ~~ ut ' > in position and 
time space seems to be rather difficult due to the various frequency dependencies involved. 
However, in a qualitative approximation, we may replace uj by 1/t. This shows that the 
pairs propagate according to 

L 2 (l + oTi t +•••)= L 2 + aD t + • • • , t < rr 1 
(/?-(/■)) - V, A{ a/t)t ~ { D [\n(T 1 t)}- d t , rr 1 < t < r c (2.39) 

D [In L sy8 ]~ d t , r c <t 

11 



We have used (see |2.37|) the notation D = (U 2 /Wb) L 2 ~ d , and a and a are numerical 
constants of the order of unity. For very small time scales t <C T^ 1 the dynamics is not 
diffusive. The corresponding estimate given in fl2.39|) for these time scales may formally be 
obtained by a naive combination of (|A7| ), ( 2.34 ) and ( |2.38| ), which is strictly speaking not 
justfied. However, a more sophisticated analysis for the classical propagator lZ(R,t) shows 
that this estimate is indeed correct. The physical interpretation of ( |2.39| ) is that the electron 
pairs are immediately (on time scales which are not resolved by the ERMH) delocalized on 
a scale L±, due to interaction induced random hoppings. Then a diffusion with diffusion 
constant D sets in. Eventually the diffusion constant is suppressed by the time dependent 
factor [ln(ri£)] _d which saturates at [lnL sys ] _d for t > r c . 

We also expect that, due to the frequency dependence of z/ eff (u;), the number of diffusing 
states depends on time as N(l/t) with N(u) defined in (|2.29|) . The time dependent pair 



size therefore behaves as L\ ln(Fit) for rY 1 <C t < r c and as LilnLgyg for r c < t. In 
recent numerical simulations 0,0 of two "interacting" kicked rotators, Borgonovi and 
Shepelyansky indeed observed the logarithmic increase of the pair size with time. 

As a further illustration, we consider the probability that an electron pair performs a 
transition from a product state \Rj > to another state \Rj > in the time t > 0. This 
probability is given by 



<Rj\e- tm \Rj> 



where 



1 



47T 5 



dE \ duo e 



-iwt 



(+) 



(±) 



(Gg) ~' W =< Rj\ 



m m v + %(G$-\E-% 



\Rj> 



E±ie-H 



(2.40) 



(2.41) 



denotes the matrix elements of the advanced or retarded Green's function. The ensemble 
average with respect to ( of the transition probability can be calculated from the functional 
F(J) (choosing a suitable source matrix J and taking the derivatives). Omitting technical 
details, we give the result of a perturbative evaluation in the diffusive regime using the a 
model ( P3ED : 



<Rj\e- int \Rj> 



dE / dui e 



-iujt 



X 



^vQ) 



\ 



[e-4 + ItovGj) + Qr^G) - fc) ) 



E-rf R + ±T v(j)) +{lT lV (j)-lu 
x(2n)- d fdqn(q,uj)e iq{R - R) . 



(2.42) 



The contribution in the last line is the "diffusion" propagator in position and frequency 
space with TZ(q,u) = [a c s q 2 — iui h> e fi{u))}~ 1 as above. Physically, it describes how the 
electron pair diffuses along its center of mass coordinate R [with the modifications shown in 
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( p.39j )1- The only information about the relative coordinate j is contained in the Lorentzian 
prefactors in ( J2.39D . These effectively define a resonance condition for the diagonal energies: 
Wr ~ V R \ ~ r lining (j),v(j)]. For the well coupled product states (v(j) ~ 1, v(J) ~ 1), 
we recover that essentially states inside a band of width r x are coupled by the diffusive 
transport processes. For the badly coupled product states (v(j) ~ e~ 4 ' Jr ' //Ll , v(j) ~ e~ 4 ' J '/ Ll , 
with \j\, \j\ ^> L\) this resonance condition is practically never fulfilled and the transition 
probability becomes exponentially small due to the Lorentzian factors. However, for certain 
improbable realizations of the rf R the resonance condition may be fulfilled even for two badly 
coupled product states. Then the corresponding transition probability experiences a strong 
enhancement. 

F. TIP Level Statistics at L > L x 



The a model action (|2.36| ) also provides a very convenient tool to discuss level correla- 



tions. Actually, the two point correlation function, defined by 

UUJ) = (p(E)) ( P \e + <j)) ( {p{E)) {p{E + U)) ~ {p{E) P{E + U)) ) ' (2 - 43) 

(p(E) is the unaveraged total density of states) can be calculated as in Ref. ||. Of course, in 
this way, we only obtain the correlations of the well coupled (diffusing) pair states. The iso- 
lated pair states are mainly uncorrelated with the well coupled states and among themselves. 
Therefore, their contribution to the correlation function (|2.43|) essentially cancels out. For a 
finite system of size L, the diffusive dynamics defines an energy scale E^\ the "pair Thouless 
energy", which is in our case determined by the implicit equation V e ^{E^)/ L 2 = EJ?\ For 
frequencies smaller than E^ the second term of the action /^[Q] i n ( |2.36|) dominates the 
level correlations, ^(w) can exactly be expressed by the same integral (over one Q-matrix) 
as in |J. Therefore we recover the random matrix result 

Y 2 {u) = Yf° E) (-^-) , AH" 1 = L d u«{u>) , M<S« (2.44) 



with v e f[ as in ([2.371) , ( |2.38|) . A(c<j) is a frequency dependent effective levelspacing and 



F 2 (r) is the universal spectral correlation function of the Gaussian orthogonal ensemble 
(r is the energy difference measured in units of the level spacing). For higher frequencies, 
i.e. E^) <^C \uj\ <ti T 1 the first (kinetic) term of /^[Q] has to be taken into account. The 
perturbative evaluation of l^(^) is completely equivalent to the diagrammatical calculation 
of Y 2 (uj) given by Alt'shuler and Shklovskii [J33] . We get 

1 / \ d/2 

u 2 \V eS (uj) j L 2 J 

Here V e g(uj)/L 2 can be viewed as a frequency dependent Thouless energy that sets the scale 
for the level correlations in the diffusive regime. 

In summary, we state that the field theoretical o model formulation of the effective 
random matrix model enabled us to recover many well known results for the localization 
problem of noninteracting particles. 
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III. TIGHTLY BOUND PARTICLES IN ONE DIMENSION 

In this section, we are concerned with a particular example in one dimension which allows 
for a more microscopic treatment than the mapping onto the effective random matrix model. 
Our aim is to show that the physical implications and results found in the last section are 
not restricted to the random matrix approach and are indeed generic for the two interacting 
particle problem. Moreover, in the present case we are able to generalize our treatment to 
more than two particles. 

The essential idea is to compare the microscopic formulation of the problem of two 
interacting particles in a Id random potential with the related problem of independent 
particles in a quasi Id disordered strip. To this end we start with the case of two interacting 
particles and consider the Hamiltonian 

H = - l -{d 2 Xi +d 2 X2 ) + V{x l ) + V{x 2 ) + U{\x l -x 2 \) , (3.1) 

where X\ and x 2 are the coordinates of the two particles, U(\x\ — x%\) is the interaction 
potential, and V(x) the Id random potential with 

(V(x)) = , 
(V(x)V(x')) = c 1 5(x-x') . (3.2) 

Introducing new coordinates x = (x\ + x 2 )/2 and y = x\ — x 2 we have 

H = -\d 2 x + l -H T {y) + V(x + y/2) + V(x - y/2) (3.3) 

with the transverse Hamiltonian 

H T {y) = -2% + 2U{\y\) . (3.4) 

For later convenience we define a rescaled Hamiltonian H = 2H. We will study the transport 
properties of H at the two-particle energy E (i.e., at E = 2E for H) and define a one-particle 
reference momentum kp by E = 2k 2 F /2. Next, we diagonalize the transverse Hamiltonian, 

H T (y)<f> n (y) = (-2d 2 y + 2U(\y\)<p n (y) = ^Mv) , (3-5) 

where the <p n {y) are real eigenfunctions and we restrict ourselves to the even case, i.e. 



0n(y) = 4>n{—y)- Similar to Dorokhov |24] we assume an attractive potential between the 
two particles. In particular we consider a bag model, where £7(|y|) = oo for \y\ > L and 
[/(M) = otherwise. For this example the transverse eigenfunction system is simply given 
by 

0n(y) = —r= cos(k n y) , 

k n = ^r(n + -) (n = 0,1,2,...) , 

e n = 2k 2 n = 2{^j\n+ 1 -) 2 . (3.6) 
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For the case of two electrons in a random potential an attractive interaction is of course 
unrealistic. However, analytical arguments [1^,1^,^] and in particular numerical studies 
based on the transfer matrix approach ||T6| , |18|| indicate that the sign of the interaction is 



basically immaterial for the transport properties of the system. In j[8] no essential difference 
between a repulsive Hubbard-type interaction and an attractive bag interaction was found. 
Therefore we feel justified to proceed with the conceptually much simpler case of an attractive 
bag. Also, one might think of a particle pair which is indeed (tightly) bound, e.g. an exciton. 
We project H on the transverse eigenfunction system, 



H„ 



V n.m. \ X I 



\h\K) 



~lfl + e n 



Onm i v nrn \X) 



dy <j>M<t> m (y) [V(x + y/2) + V{x - y/2)} . 



(3.7) 



Thus we have mapped the problem on a Id, many channel Hamiltonian. We define a diagonal 
matrix of longitudinal momenta, 



k = diag(«4 . . . K r 



K T 



2(E - e r 



(3.8) 



The N open channels of the problem are defined by the condition k n = Jt n /2 < kp. 



A. Mapping onto a Quasi-ld Transfer Matrix Approach Without Interaction 

For comparison, let us now consider one electron in a disordered 2d strip, 

1 



H 



M 1 +d 2 J + U(\x 2 \)+V(x 1 ,x 2 ) . 



(3.9) 



Here, X\ and x 2 are the two space coordinates, t/(|x 2 |) is the confining potential defining 
the strip and V(xi,x 2 ) the 2d random potential with the statistical properties 



(V{ Xl ,x 2 )) = , 

{V(xi, x 2 )V(x' 1 , x' 2 )) = c 2 5(xi - x[)S(x 2 - x' 2 ) 



(3.10) 



For notational simplicity we do not explicitly distinguish between the two cases considered 
here and hope that no confusion will arise. In close analogy to the previous problem of two 
particles we diagonalize the transverse Hamiltonian 



H T (x 2 ) = --% + U(\x 2 



(3.11) 



so that 



H T (x 2 )(f) n (x 2 ) = e n (j) n (x 2 ) 



(3.12) 



Again, we restrict ourselves to the even case <p n {x 2 ) = <fi n {—x 2 ). If we choose a box of width 
2L with infinite walls as confining potential the transverse eigenfunction system is given by 



!6|). Projection of ( p.9|) on the transverse modes gives 
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H„ 



Vn 



■ / n\- 1 - 1 \Ym 



) 



-~^ + e r 



Onm ~r Vnm\%l) 



dx 2 (j) n {x 2 )(j) m {x 2 )V(x X ,X 2 ) 



(3.13) 



A matrix of longitudinal momenta can be defined as in ( ft.8| ). The reference momentum kp 
of the two electron problem can now be interpreted as the one electron Fermi momentum 
and the open channels are characterized by the condition k n = v / 2e n < k F . Comparing ( p. 71) 
with (|3.13|) we see that both the strip problem and the two electron problem are described by 
very similar quasi Id (many channel) Hamiltonians. The only important difference resides 
in the definition of V nm (x) in the two cases. In the next step we apply the transfer matrix 



method to (|3.7|) and (|3.13| ). 

We define a longitudinal iV-channel wavefunction ip(x) T = (ipi(x), . . . ,iPn(x)) so that 
the Schrodinger equation for a Hamiltonian of the type ( |3.7j ) or ( |3.13| ) reads 



Eif){x) = Hip(x) 



■-^"{x) + eip(x) + Vip(x) 



(3.14) 



where e is a diagonal matrix with (e) nn = e n and the primes in ip"(x) indicate the derivative 
with respect to x. We introduce a transfer matrix A by 



ip(x) 
ip'(x) 



1 

2(e-E + V) 



ip(x) 
if;'(x) 



= A 



ip(x) 
ijj'(x) 



(3.15) 



Redefining the wavefunction through 



v(x) 



%1)\X) + IKlfjiyX) 



C 



ip(x) 



(3.16) 



y/2n \ ^ / ( x ) — iKijj(x) J \ ip'(x) 

[we recall that k is the diagonal matrix ( |3.8|) ) we have v'(x) = CAC~ 1; u(x) = T(x)v(x) with 

( k \ if V(x) -V(x) \ 1 



T(x) 



\0 -k) ^\V{x) -V(x) j 
To get rid of the first matrix on the r.h.s. of ( |3.17| ) we redefine the wavefunction again, 

~, x f e lKX \ , , 



(3.17) 



v\x) 



e 



v\x) 



and finally arrive at 



v'(x) = T(x)v(x) 



with 



T(x) 



e lKX V{x)e l 
e iKX V(x)e iK 



c V{x)e- 
■ ; V(x)e~ 



(3.18) 



(3.19) 



(3.20) 



It is not difficult to see that ( |3.19| ) is equivalent to the following evolution law for the full 
transfer matrix M(x), 
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M(x + 6x) 



JM 



M(x) 



(3.21) 



where the "building block" 6M is given by 



5M 



x+Sx 



T(x')dx' 



8v\ —5v 2 
5v\ —8v\ 




(3.22) 



with 



8v\ 



5v 2 



x+Sx 



dx'e~ lKX V(x')e lKX 



X + &X 



dx'e~ lKX V{x')e~ tKX 



(3.23) 



Under the crucial assumptions that both the blocks of T(x) in (|3.20|) and all the ma- 
trix elements within these blocks are statistically independent, we deal with the so-called 
isotropic situation and the evolution law ( |3.21|) leads to the Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equation [IT] for the eigenvalue distribution of M(x). This equation accurately 
describes the transport of independent electrons in quasi Id wires and is — at least in the 
unitary case -- completely under control [|T^,|TB[]. In the case of a 2d strip of width not 
significantly exceeding the elastic mean free path l e \ the exponential factors e ±lKX in (|3.20| ) 



oscillate rapidly so that the assumption of statistically independent matrix elements seems 
very reasonable. This argument provides a microscopic justification for the application of 
the DMPK equation to quasi Id wires [11]]. A main purpose of our present study is to point 



out that the DMPK equation also describes the transport of two interacting electrons in a 
disordered Id system, provided the pair size does not exceed l e \. This statement should be 
practically obvious since both the strip and the two electron problem have been formulated 
in almost identical technical terms: The DMPK equation is as well-justified for the latter 
problem as for the former. 

B. TIP-Localization Length 



Once the applicability of DMPK theory is guaranteed the next step is to identify the 
localization length of the problem in terms of the given system parameters. The connection 
between the building block Q3.22J ) and the localization length £ in the DMPK approach is 
given by 



N 2 

—5x = (tr(¥b)) 



(3.24) 



i.e. 



—rox 



x+5x 



dx\dx 2 (tr 



—e 
Ik, 



-in(x\—X2) 



V( Xl )-e 



_ -in(x 1 -x 2 ) 



H 



V(x 2 ) 



(3.25) 



Any further evaluation of (|3.25|) depends on the precise form of V(x) and there, finally, the 
differences between the two cases studied here come in. 
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Let us begin with the disordered strip. Inserting V(x) as defined in (|3.13| ) in (|3.25| ) we 



get 



—^ox 



x+&x rL 

dx\dx 2 \ dy\dy 2 tr 

J-L 



1 

— ( 

K 



-e- tK{xi - X2) ^y 1 )cl)(y 1 ) T -e- l ^- X2 ^(y 2 ) ( j>(y 2 " T 



(V(x 1 ,y 1 )V(x 2 ,y 2 )) 

and after using ( |3.10| ) and then ( |3.6| ) 

N 2 r L 1 

Sx = c 2 5x / dy^2 (Pn(y) 2 (pm(yf 

J—L nrri Hn^rn 

(-2 



(3.26) 



i 



nm lv W^m 



^* J nm K n Kr,„ 

Evaluating the remaining sums approximately (see appendix 



(3.27) 



N-l i 

n=0 ^n 



L 

~2 



N-l 



2-~i K 2 

n=0 ^n 



1 L 2 \n2N 



2vr 2 N 



(3.28) 



we can neglect the diagonal contribution in ( |3.27| ) as a 1/N correction. Together with (see 
appendix 0) 



TT^ 1 ki~ d 



Cd 



for d = 2 we finally arrive at 

N 2 n k F L 



d kpl e i 



& £ = —Nl el 

1T Z 



(3.29) 



(3.30) 



e 16 l d 

This is the well-known relation between the localization length and the elastic mean free 
path in quasi Id disordered strips. In the following we will show that a similar relation also 
holds in the case of two electrons with a bag interaction, although the detailed reasoning is 
a little more involved than in the previous case. From ( |3.7|) and ([3.25 ) we have (instead of 

N 2 



£ 



-Sx = 4 



x+Sx 



dx\dx 2 \ dy\dy 2 tr 



-iK(x\—X2) 



K 



$(vi)<f>(yi) - e 

K 



-ik(xi—X2) 



0(^2)0(1/2 



{[V(x l + 3/1/2) + V( Xl - yi /2)] [V(x 2 + y 2 /2) + V(x 2 - y 2 /2)\) 



(3.31) 



According to (|3.2|) the average over the random potential leads to four 5-functions so that 
x\ — x 2 = (s±yi + s 2 y 2 )/2 with s\ j2 = ±1 (independently). Since the typical scale for 5x is 
l e i the ^-function condition can be always fulfilled and we arrive at 

72 



N r / 1 r 

—5x = A Cl Sx J2 tr (- d Vl e- iKSl ^ l/2 0( ? / 1 )0( ? / 1 ) T e- iKSl ^ /2 

~ J dy 2 e-^y^^y^iy^e-^^ 2 



(3.32) 



or 



2 



—— = 16ci tr 



l -B l -B 

K K 



(3.33) 



where 



B nm = I dy cos((K n + K m )y/2)<f>n(y)<f>m(y) ■ (3.34) 



L 



Employing ( p.29|) for d = 1 and introducing r\ n = J N 2 — (n + 1/2) 2 , Q3.33J ) can be rewritten 



as 



rr^^r 131 - ■ (3 - 35) 



In appendix |D| it is shown that (|3.35|) leads to the final result 



^WfRiTTl) ' (3 ' 36) 

Using N = kpL/ir and L p» l e \ (which is the maximal pair size allowed within the present 
technical framework) we have 

e~M e 2 i , (3.37) 



a result which has first been derived by Dorokhov p4j for two particles interacting via 
a harmonic potential. Eq. ( |3.37| ) shows that the localization length for two interacting 
particles is significantly enhanced over the one-particle value £ ~ l e [. 

C. Scaling of the TIP-Conductance 

In conclusion, we have demonstrated the validity of the DMPK equation for the transport 
properties of two electrons in a narrow bag in Id. As our main argument we have pointed 
out the almost complete analogy to the the case of a disordered quasi Id strip. Identifying 
the localization length in terms of the system parameters we recovered Dorokhov's result 

As probably the most important consequence of our investigation we can identify, as in 
the a model formulation above, a scaling parameter in the two electron problem, namely the 
conductance g = 2tr[(M^M + (M^M)^ 1 + 2) _1 ] of the associated strip problem. Although 
this quantity does not seem to have an immediate physical interpretation for two interacting 
electrons and should perhaps be regarded as some sort of auxiliary variable it is nevertheless 
of considerable value: Its length dependence in the whole range from the diffusive to the 
localized regime is under control |I(J and the crossover point between these regimes is given 



by g(0 = 1. In this way we have again justified the assumptions that underly Imry's 
application [|I7]] of the Thouless block picture to the problem of interacting particles in a 
random potential. We think that Imry's pair conductance gi should be identified with our 
conductance g of the associated strip problem. This identification explains the remarkable 
success of the scaling block picture in the present context and clarifies the nature of the 
scaling parameter. 
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D. Generalization to Many Particles 

The above considerations for two interacting particles in a random potential can be 
generalized to larger particle numbers provided the particle density (in the case of Fermions) 
is low enough to neglect effects of the Pauli principle. In this paper we refrain from discussing 
the full many body problem, where quasiparticles are the relevant degrees of freedom. A 
particularly important issue in this context would be the lifetime of interacting quasiparticles 
as a function of their excitation energy. For a first numerical study of the influence of the 
Pauli principle on interaction-assisted coherent transport see [[34]. 

Let us consider the M-particle Hamiltonian 

1 M _ 

H = -^Y,{dl + V(x l ))+Y,U(\x l -x J \) , (3.38) 

i=X i<j 

where U(\x — x'\), as before, represents an attractive bag of size 2L. Introducing the center 
of mass coordinate x and M — 1 relative coordinates, 

x M 

2 = 1 

Vj = x j ~ x i+i (; = 1, . . . , M - 1) , (3.39) 

we get 

i M 

H = -TTyr d l ~ f( d Vn • • • : d VM-l) + Yl V ( X + 9i{Vl, ■■, VM-l)) 

i=X 

+ Y, U ( h i j (.Vi,...,VM-i)) • (3-40) 

i<j 

Here, /, gi, and h%j are functions which we do not have to specify in detail for our subsequent 
considerations. In complete analogy to the previous section we can define a transverse 
Hamiltonian and its eigenfunction system by 



H T = ~f(d yi , ■ ■ ■ ,d VM _i) +J2 u ( h ij(yi> ■ ■ ->Vm-i)) , 

i<j 

H T ip ni ... nM _ x {y u . . . ,Vm-i) = e ni ... nM _ 1 (f ni .„ nM _ 1 {y u . . .,y M ~x) ■ (3.41) 

We do not specify any particular symmetry for the transverse eigenf unctions 
( Pni...nM-i(yij ■ ■ ■ i Vm-x)- Projecting the full Hamiltonian H on the transverse eigenfunc- 
tion system we arrive at 

r m...n A /-i I o2 i , \ r r i t>"i— «m-i , 



2M 
where 



hzZ:: = -^ + wJ^ . . . w.^ + y^Z::^ . ( 3 - 42 ) 



V n'!.'.'ZJ( x ) = J dyi... dy M -x <p ni ... nu -i(yi, ■■■, Vm-x) 

( S V ( X + 9i{VU ■ ■ ■ > VM-x)) J ^n' 1 ...n^_ 1 (VX, ■■■, VU-A 



(3.43) 
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Defining the "building block" SM as in ( |3.22| ) and keeping in mind the relation ( |3.24|) we 
have 

(ti(tfb))= I* * dxxdxi hr(-e- l ^~ x ^V{ Xl )-e- iK{x '- X2) V{x 2 ))\ , (3.44) 



where K ni _„ nM _ x = J2(E — e ni ... nM _ 1 ). The relevant average in ( p.44[) is straightforwardly 
calculated, 

/ M M \ 

( X v ( x ± +9i(vi--- vm-x)) Yl v ( x ? + 9j(y'i ■ ■ ■ v'm-i)) ) = 

k 3 M 

2jfrJ2 s ( x ^- x ^ + 9i(yi---VM-i)-9j(y' 1 ---y' M -i)) • ( 3 - 45 ) 

/ e( i t j=l 

Both the %\ and the x 2 integration in Q3.44J ) can be easily performed and we get 

k 3 M r 

(tr(6 f 6)) = F 8x ^2 d y 1 ... dyM-idy'x ■ ■ ■ dy' M _ x 

tr f- e - iK 0n-Bi)(p(y u . . . , yu-iMvi, • • • , Vm-i) T 



e-^-^ V (y[,...,y' M ^My' 1 ,...,y' M _ 1 ) T 

„3 



-p-8xtr (-B-bA (3.46) 

2kpl e l \K K J 



with 



<Pn 1 ...n M _ 1 {yi,---,yM-l)<Pn' v ..n' M _ 1 {yi,---,yM-l) , (3.47) 

which corresponds to (|3.34|) . For further progress we have to rely on certain estimates. 
Let us assume that the phase factor exp(i(« ni ... nM _ 1 + K n > ...„' )9i(yi, ■ ■ ■ , Vm-i)) fluctuates 
strongly when the yi vary between — L and L. Then it can be replaced by its average denoted 
by c p h and we have 

Our assumption is only correct if «m...n M _ 1 ^ 1/L. This means that channels with high 
transverse excitation energies and small longitudinal momenta K ni ...n M _ 1 are not correctly 
described by our approximation. Comparing the result of our estimate with the more ac- 
curate calculation for M = 2 in the previous section we can determine the error incurred. 



(3.49) 



3.24J) and (|3.48) we have 




1 r 2 , k 3 

£ " N* 2k F l el ni ^ M _ 


1 

.1 K m...n M -i 
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We estimate the sum in (|3.49|) by an integral approximation, 

£ ^— K f *' (£) 2(ib) • (350) 

m...n M -i ^m.-.TiM-i J0 z ^ fc J 

where E is the M-particle energy at which we study the transport problem, and p(e) is the 
density of transverse energy levels. To continue, we replace p(e) by its average p and assume 
that E — e ~ A (with A the transverse level spacing) provides the cutoff for the otherwise 
divergent integral in ( |3.50| ): 



E 



1 /•! p(x) 



de p(e)—— — - dx 



2(E-e) 2 Jo 1-x 
p z -1-5 dx 



2 Jo 1-x 
= -^ln<5 . (3.51) 

With E ~ kp, the number of transverse channels N = (kpL) 1 , p = N/E, and 5 = A/E = 
1/iV we finally arrive at (disregarding numerical prefactors) 

1 (k F L) M ~%, 
£ ~5ln[(*,L)«-] ■ (3 - K) 



For the case M = 2 we almost recover the result ( |3.37|) of the previous section. The 



logarithmic suppression in fl3.52| ) is an artifact of our approximations. It arises precisely 



from the levels with E rj e, i.e. k ~ 0, which are not well-described by our approach. 
Therefore we discard the In-correction in (|3.52|) and have as our final result (for L w l e i) 



Z~(k F l el ) M -H el . (3.53) 

This shows how the enhancement factor for coherent multiple particle propagation increases 
with particle number M. 



IV. CONCLUSION AND DISCUSSION 

In this work, we have developed and further improved two analytical approaches to the 
phenomenon of interaction-assisted coherent pair propagation. One, the a model formula- 



tion based on an effective random Hamiltonian, was inspired by Shepelyansky's view |L6 
The other, a transfer matrix approach to two or many tightly bound particles, was built 
on the pioneering ideas of Dorokhov |24| . In its most succinct form, the message from our 
investigations is the following: Transport of interacting pairs or of small aggregates of parti- 
cles can be described in terms which have been introduced in the context of, and are therefore 
familiar from, one-electron theory. 

In particular, pair transport is governed by a scaling parameter, the effective two-particle 
conductivity. This is a rigorous statement, once the models introduced above are accepted, 



and adds to our confidence in Imry's generalization [I7| of the Thouless scaling picture. 
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Furthermore, there is an energy scale, the two-particle Thouless energy, associated with dif- 
fusive pair transport. It corresponds, as in one-electron theory, to the time it takes to diffuse 
through the sample and defines the scale, for which deviations from random matrix theory 
occur in the spectral statistic of diffusing pair states. The precise pair dynamics, however, 
differs from the one-particle case in that both the diffusion constant and the density of states 
are frequency dependent. This is a clear manifestation of the pairs' composite character, 
since it is connected with a logarithmic growth of the pair size and a corresponding loga- 
rithmic attenuation of the diffusion constant as a function of time. Apart from confirming 
important aspects of these conclusions, our second (transfer matrix) approach enabled us 
to study, in a rather general form, the behavior of more than two tightly bound particles. 
It turned out that the enhancement factor, i.e. the ratio between the M-particle and the 
one-particle localization length, grows as the (M — l)th power of the latter. 

Of course, both models introduced in this paper have their shortcomings. Foremost, our 
a model formulation relies heavily on the assumption that the effective Hamiltonian catches 
the essential physics of the problem. We can imagine two types of criticism: (i) quantitative 
statistical assumptions like the size of the fluctuating interaction matrix elements or the 
form of the decay towards the band edges are incorrect or unrealistic. Corresponding modi- 
fications would lead to quantitative corrections without precluding our whole approach, (ii) 
Qualitative assumptions like the statistical independence of the matrix elements are incor- 
rect in an essential way. This would call into question the a model formulation as such, 
since it is hard to imagine how a similar a model can be obtained in the presence of impor- 
tant correlations. One should also note that the one-particle localization length L\ is the 
smallest length scale resolved by the effective Hamiltonian. For length scales L < L\ effects 



associated with the transport of single particles have to be taken into account [pq| . In this 
case the interpretation of the two-particle conductance is less straightforward. 

The major drawback of the transfer matrix model, on the other hand, is its restriction 
to tightly bound particles. This restriction is necessary to prevent the disorder from being 
relevant for the relative motion of the particles. Otherwise relative and center of mass motion 
would have been coupled and the problem much harder to solve. Naturally, however, the 
requirement of strong binding restricts the physical situations to which the model may be 
applied. 

Combined, the two approaches explained in this paper represent our analytical knowledge 
concerning interaction-assissted coherent transport. For the future, we believe that it is an 
important and worthwhile task to devise and investigate an effective Hamiltonian for more 
than two particles. This would be another step towards understanding many-particle effects 
and it would be free of the constraints inevitably associated with the transfer matrix method. 
It is also of interest to investigate the question, whether the effect discussed here is related 



to a mechanism for the enhancement of persistent currents proposed some time ago |36| 
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APPENDIX A: EXPANSION OF THE EFFECTIVE POTENTIAL ( g]3g ) 



In this Appendix, we derive an expression for the effective Q-potential Jr(Q) which is 
more suitable for the comparison with the standard a model of a disordered metal C met [Q}. 
For this, we consider a quantity like 

/ = str ln(a + bk + cQ) (Al) 

with some complex parameters a, b, c and A, Q as in ( |2.32| ). In the grading for advanced 
and retarded Green's functions, A takes the form A = L J. The quantity / only depends 
on the radial parameters 9 || of Q = T~ 1 kT, i.e. we may choose T = exp(8ai), where 
o\ = L j and 9 is 4 x 4 supermatrix with two or three independent radial parameters || 
(depending on the symmetry class). Then we have a\ka\ = —A and <j\Qo\ = —Q. Inserting 
1 = o\ and permuting the matrix products in the argument of the logarithm in (JAl|), we 
thus obtain: 

/ = str ln(a -bk- cQ) = -str ln((a + bk + cQ)(a -bk- cQ) 



a 2 - (6 + c) 2 -be Qk + kQ-2 . (A2) 



1 
= -str In 
2 

We now expand J, for the case in which either \bc\ -C \a 2 — (b + c) 2 | is valid or the matrix 
AQ = (Qk + kQ — 2) may be considered as small. The latter condition is clearly valid for 
a perturbative treatment where Q can be parameterized H as Q = A"|~^P — A(l + 2iP + 

(2iP) 2 H ) with APA = -P. Then AQ = 8(iP) 2 + C(P 4 ) is indeed a small parameter 

and expanding (|A2|), we find: 

be 

str(QA) . (A3) 



a 2 - (b + c) 2 
Applying (|A3|) on fn(Q) defined in (|2.32|) , we obtain: 




Here, the terms with \j \ < L c ~ L x In L sys can be considered as selfaveraging with respect to 
the random variables rj^, due to overlapping resonances. The typical (most probable) values 
of the remaing terms (with \j\ > L c ) are exponentially small since v(j) ~ e~ 4 ^'' Ll . They 
can be neglected whereas the other terms may be replaced by their //-average. In addition, 
we assume that the integration range 2W& for r\ is much larger than the other energy scales 
u>, r , r x and that the //-density po( 7 ?) varies rather slowly on these "small" scales. The 
77-integral may then be extended to all rj with — 00 < 77 < 00 and po(rj) can be replaced by 
Po(E). Performing the ^-integration, we find: 
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f R (Q) ~ -ijco h (j±) Po (E) str(QA) 



with 






•£i _ u(j) 

|j]<L c W w 



(A5) 



(A6) 



In order to obtain a more explicit expression for this function, we would have to specify the 
choice for v(j). However, the relevant limiting cases are: 



10 ) on V on 



. s- = B"0)r . r 1 «|^| 






4 lb, 



r. 



-i < 



w|<Ti , 



/i 



a; 



£i 



-,d 



In (rir c 



M < T" 



(A7) 

(A8) 
(A9) 



Here, r" 1 = A cff = Fiv(L c ) is the effective level spacing of the well coupled pair states, 
introduced subsequent to Eq. ( |2.28| ). The cutoff in the sum (|A6|) is either given by \j\ < 
L eS (oo) [cp. (|2.29| )1 or \j\ < L c , resulting in (|A8|) and (|A9|) , respectively. 



APPENDIX B: SOME ESTIMATES 

In this appendix we estimate the discrete sums X) n (l/« n )> Sn(l/ K n) contributing in 
Q3.27I ). With IE « N 2 tc 2 /L 2 , k 2 n = n 2 (n + l/2) 2 /L 2 « n 2 n 2 /L 2 , and approximating the sum 
by an integral we have 



Analogously, we get 



JV-l 1 

E- 

71=0 K ™ 



JV 


L 



dn 



2E - hi 



N 



dn 



ttN Jo ^1 _ n 2/N 2 
L r 1 dx 



vr Jo v^T" 
L/2 . 



r/ 



N-l 

E 



L 2 



n=0 K n K 2 N Jo 1 



1 dx 



Xr 



(Bl) 



(B2) 



which is logarithmically divergent. However, we have to take into account the energy dif- 
ference between E and the transverse energy of the last channel. This energy gap is of the 
order of the level spacing and hence of relative magnitude 1/N. Therefore, 
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^ 1 L 2 /"i-i/iv dx 

L 2 1, 1 + x l-i/JV 



:-ln- 







2tt 2 AT 
This completes the derivation of ( |3.28| ). 



w 2 N2 1-x 

L 2 ln2iV , s 

(B3) 



APPENDIX C: NORMALIZATION OF THE RANDOM POTENTIAL 



To estimate the constants c\ and c 2 appearing in ( |3.2| ) and ( |3.10| ), respectively, we recall 
the normalization condition for models with random white noise disorder ||, 

(V(x)V(x>) = -^—5(x-x>) , (CI) 

Ztiut 

where v is the local density of states and r the elastic mean free time (we put H = m — 1). 
With r = lei/kp and iV = (kpL/TrY in d dimensions we get for the total density of states 

>=% = $**& («) 

and hence 



ciE " 2 f \nj 



so that 



„ = A = ^-^ (C3) 



I _d-l 1,4-d 

27TZ/T O /CF'eZ 



which is exactly ( |3.29| ). 



APPENDIX D: FINAL EVALUATION OF 1/f 

We have to estimate the r.h.s. of ( |3.35| ). Using ( |3.6| ) and ( |3.34| ) an elementary calculation 
shows that 

1 „ sin \{K n + K m )L/2 + gi7r(n + 1/2) + s 2 ir(m + 1/2)] 
nm ~ 2 Sl ^ ±1 («n + «™)L/2 + si7r(n + 1/2) + s 2 n(m + 1/2) ' [ j 

With the abbreviation f(x) = sin(x)/a; ( |3.35| ) can therefore be rewritten as 
\ = Y EE — / (^ [^ + ^ + s i(^ + V2) + s 2 {m + 1/2)]) 

f (^[Vn + Vm + s 3 {n + 1/2) + s 4 (m + 1/2)]) . (D2) 
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With the help of the approximation scheme 

N - X 1 rN dx 



£ ±-h{n + 1/2) « J 



VN 2 - x 2 
1 ds 



h(x) 



h(Ns) 



o VT 

= [ n dip h(N simp) (D3) 

Jo 

for an arbitrary function h(n + 1/2) we find 

1 1 f w / 2 

- = — 22 d<pid<p 2 /(ttN [cos if 1 +cos^ 2 + sisin^i + s 2 sin^2]) 

/ (7rA^ [cos^i + cos(y9 2 + s 3 sin^x + s 4 sin<^2]) • (D4) 

The main contributions to the integral arise from the regions where [. . .] ~ 1/N. Also, 
the two factors in the integrand should be in phase, i.e. s\ = S3 and s 2 = S4 (diagonal 
approximation). Therefore we are left with the expression (using sin p> + cos ip = sin(ip + 

t/4)A/2) 

1 1 / >7r / 2 / TlN \ 

~i = ~ki J-n/2 difld(P2 ^ I V2 [Sil1 ^ 1 + 7r/4) + Sil1 ^ 2 + ^ J ' (D5) 

for which the integration range can be restricted to a narrow strip around the line <pi = 
— 7r/2 — ip 2 . A first order Taylor expansion in the deviation Aip from this line yields 

1 1 rO poo /ttN \ 

7«r dtp d(Ap)f 2 -=cos(^ + 7r/4)A^ . (D6) 

£ t e j J-k/2 J-infty \y/2 J 

With the help of the integral 

00 , sin 2 (ax) n 



dx . \/ = — (D7) 

axV a 



we get 



7T 







aV 



v/2 



£ Z e / J-k/2 ' TTN J COs(<y9 + 7r/4) | 

\/2 W 4 aV 



M eZ tan(7r/8) M ei 
which is exactly (|3.36|) . 
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iV7 e z J-tt/4 COS 99 

v^ ln tan(3V8)^V2 M1 + ^ ) _ 
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